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FINAL  REPORT  ON  IN-HOUSE  EXPLORATORY 
DEVELOPMENT  PROGRAM  ON  MICROSTRIP  ANTENNAS 


INTRODUCTION 

Due  to  its  adaptability  to  a wide  variety  of  aircraft,  missile, 
projectile  and  phased  array  applications,  the  microstrip  patch 
antenna  has  become  a subject  of  great  interest  to  antenna  and  systems 
designers  through  the  scientific  community. ^ It  is  of  particular 
interest  in  areas  where  the  thickness,  cost,  and  reliability  of  the 
antenna  are  of  prime  importance  such  as  in  projectile  fuzing  applica- 
tions. Stripline  slot  antennas,  although  equally  adaptable  to  many 
systems,  have  proven  to  be  too  costly  and  unreliable  in  many  instances. 
Much  has  been  said  and  written  regarding  the  various  forms  of  the 
microstrip  patch  antenna  and  their  potential  uses.  Reference  2 in 
particular  gives  a good  overall  view,  and  Reference  3 provides  addi- 
tional data;  however,  there  is  very  little  information  in  the  litera- 
ture which  could  be  classified  as  basic  design  data.  The  Naval  Surface 
Weapons  Center  has  undertaken  an  in-house  exploratory  effort  directed 
towards  producing  and  disseminating  a practical  design  procedure. 

BASIC  MICROSTRIP  PATCH  RADIATOR 

Figures  1 and  2 show  a microstrip  patch  radiator  etched  on  a 
teflon  fiberglass  substrate.  In  this  particular  configuration,  as 
shown  in  Figure  3,  energy  is  fed  to  the  patch  through  the  substrate. 

The  patch  can  also  be  fed  on  its  edge  as  shown  in  Figure  4,  or  notch 
fed  as  in  Figure  5.  The  latter  two  feed  methods  are  in  fact  the  most 
commonly  used.  They  allow  the  entire  feed  network  to  be  etched  on  the 
same  substrate  as  the  patch,  a definite  advantage  in  multi-patch 
applications,  an  example  of  which  is  shown  in  Figure  6. 


1.  Black,  L.  M. , and  McCorkle,  J.  W. , "Preliminary  Report  on  the 
In-house  Exploratory  Development  Program  on  Microstrip  Antennas," 
NSWC/WOL  TR  75-200,  December  1975. 

2.  Khulman,  E.  A.,  Microstrip  Antenna  Study  for  Pioneer  Saturn/Uranus 
Atmosphere  Entry  Probe,  NASA  Report  CF-137513. 

3.  Munson,  R.  E. , Conformal  Microstrip  Antennas  and  Microstrip  Phased 
Array,  IEEE  Transactions  on  Antennas  and  Propagation,  January  1974. 
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FIG.  3 DETAIL  OF  PATCH  RADIATOR  REAR  FEED  STRUCTURE 
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HG.  4 EDGE  FED  MICROSTRIP  PATCH  RADIATOR 
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FIG.  6 MULTIPATCH  MICROSTRIP  TELEMETRY  ANTENNA 
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DESIGN  PROBLEM 

There  are  two  basic  design  problems  associated  with  the  design 
and  construction  of  a microstrip  patch  antenna.  These  are  the  deter- 
mination of  the  dimensions  of  the  patch  itself,  and  the  determination 
of  the  input  impedance  at  a given  point  on  the  patch.  They  are  a 
function  of  the  thickness  of  the  substrate,  the  dielectric  constant 
of  the  substrate  material,  the  operating  frequency,  the  patch  feed 
itself,  and  ground  plane  size.  Figure  7 is  a diagram  of  the  micro- 
strip patch  and  all  of  its  parameters. 

GENERAL  THEORY  OF  OPERATION 

There  are  three  (3)  independent  modes  of  oscillation  which  can 
occur  on  the  back  fed  microstrip  patch  antenna.  A wave  can  oscillate 
between  the  ends  of  the  patch  along  the  'A'  dimension.  When  'A'  is 
close  to  a half  wave  length,  this  is  analogous  to  a simple  half  wave 
dipole  antenna. 

A second  mode  of  oscillation  can  be  set  up  along  the  'B'  dimension. 
Also,  a third  mode  can  be  set  up  along  the  pin  used  to  feed  the  patch. 
This  mode  is  referred  to  as  the  monopole  mode  since  the  pin  can  be 
thought  of  as  a top  loaded  monopole  over  a ground  plane  where  the 
patch  is  the  capacitive  top  hat. 

The  input  impedance  of  the  microstrip  patch  is  a function  of  Y 

and  all  three  modes.  Just  like  a half  wave  dipole,  the  impedance  §t 

the  center  (Y^  = 0)  is  low,  and  the  impedance  at  the  end  (Y  = A/2) 

is  very  high.  Either  or  both  modes  A and  B can  be  utilized°by  properly 

choosing  Y and  Y,  . 

, a b 

i 

i INITIAL  STUDY 

1 Initially,  the  computer  was  programmed  with  Kaloi's  equations  to 

i determine  the  input  impedance,  the  antenna  gain  pattern,  the  resonant 

! frequency,  and  the  bandwidth.  It  was  in  this  period  that  several 

corrections  were  made  to  Kaloi's  paper  (Reference  4). 

I 

Kaloi's  equations  were  based  on  there  being  only  one  mode  excited 
I to  simplify  the  calculations.  His  equation  for  determining  the 

i length  of  the  antenna  (1)  has  been  found  to  be  very  accurate. 

[ From  Reference  3: 

I 

''■(tfMi)  - + .61  (E^-11(B/H)-1155  j 


4.  Kaloi,  C.  M. , Asymmetrically  Fed  Electric  Microstriy  Dipole, 
Naval  Missile  Center,  Point  Mugu  Technical  Publication, 
TP-75-03. 
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A = PATCH  LENGTH  INCHES 
B = PATCH  WIDTH  INCHES 
H SUBSTRATE  THICKNESS  INCHES 
- FEED  POINT  INCHES 
E ° = SUBSTRATE  DIELECTRIC  CONSTANT 
7 » SUBSTRATE  LOSS  TANGENT 
AG  > GROUND  PLANE  LENGTH 
BG  = GROUND  PLANE  WIDTH 


FIG.  7 ANTENNA  COORDINATE  SYSTEM  WITH  PARAMETERS  DEFINED  AS  FOLLOWS: 
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Unfortunately,  his  equations  for  determining  the  input  impedance 
yielded  VSWR's  as  high  as  3. 5:1.0.  However,  the  sin2  distribution 
was  found  to  match  our  data  very  well.  We  combined  these  results 
into  a semi-graphical  solution  included  in  Appendix  1. 

The  ground  plane  size  was  progressively  reduced  on  several 
antennas  with  the  following  results:  (1)  As  AG  was  reduced,  the 
input  impedance  increased;  (2)  As  BG  was  reduced,  the  input  impedance 
decreased;  (3)  The  frequency  was  almost  independent  of  the  ground 
plane  size.  All  of  these  changes  were  small  enough  to  be  neglected 
in  the  final  design  equations. 

Another  trend,  which  was  small  enough  to  be  ignored  in  the  final 
design  equations,  was  an  increase  in  the  resonant  frequency  as  Yq 
increased . 

COMPUTER-AIDED  STUDY 

Having  already  acquired  data  on  a dozen  antennas,  the  computer 
was  enlisted  to  keep  track  of  the  data.  Curve  fitting  was  used  to 
augment  Kaloi's  equations  such  that  reliable  input  impedances  could 
be  predicted. 

A large  data  base  is  required  to  get  a reliable  curve  fit,  so  we 
started  building  antennas.  Eventually,  antennas  were  made  and  data 
was  taken  on  four  frequencies  from  400  MHz  to  2GHz,  on  boards  that 
were  1/8,  1/16,  and  1/32  inches  thick,  on  boards  made  by  three 
different  manufacturers,  and  with  A/B  ratios  of  .67  to  2.0.  In  all, 
over  60  antennas  were  built  and  tested  with  each  succeeding  set  of 
antennas  based  on  results  obtained  with  all  of  the  preceding 
antennas . 

Ten  (10)  programs  and  eleven  (11)  subroutines  were  eventually 
written  and  used  in  the  process  of  the  study.  As  a result  of  the 
subroutine  oriented  organization  of  the  programs,  other  users  can 
easily  calculate  any  antenna  parameter  required.  A complete  listing 
of  all  of  the  progreims  is  included  in  Appendix  2. 

FINAL  RESULTS 

Based  on  approximately  60  antennas,  the  maximum  VSWR  was  1. 3:1.0 
and  the  average  VSWR  was  1.1: 1.0.  Program  MICROAN  allows  anyone  to 
design  a complete  back  fed  microstrip  antenna  via  a teletype  terminal. 
An  example  run  is  shown  in  Figure  8 and  a flow  diagram  for  MICROAN  is 
shown  in  Figure  9.  The  only  limitations  to  this  program  are  that  the 
conditions  H < 1/8  and  A/B  > 1 be  met. 

ADDITIONAL  WORK 


Antenna  patterns  were  measured  on  several  antennas.  Very  good 
agreement  between  the  calculated  patterns  and  measured  patterns  was 
obtained . 
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Based  on  nine  notch-fed  antennas,  a graph  has  been  made  which 
can  be  used  to  augment  the  sin^  distribution  of  the  antenna  input 
impedance.  Upon  further  development,  this  could  be  incorporated  into 
program  MICROAN  so  that  notch-fed  antennas  could  be  designed  via  the 
teletype  terminal. 
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APPENDIX  A 
Size  Calculation 


C.  M.  Kaloi's  equation  for  determining  the  length  of  the  antenna 
has  been  found  to  be  very  accurate.  The  basic  problem  in  finding 
a concise  deisgn  procedure  for  determining  the  patch  dimensions  is 
the  fact  that  'A'  is  a function  of  'B'  as  seen  in  Equation  (1)  below. 

From  Reference  3: 


As  a result,  a design  procedure  has  been  set  up  whereby  the 
ratio  of  'A'  to  'B'  is  arbitrarily  chosen,  followed  by  the  calcula- 
tion of  both  'A'  and  'B'.  The  derivation  of  the  procedure  is 
explained  on  a step-by-step  basis  below. 

The  design  procedure  is  based  on  frequency  scaling  and  on  scaling 
with  the  square  root  of  the  dielectric  constant.  The  frequency  scaling 
procedure  results  from  an  approximation  in  the  first  term  of 
Equation  (1)  as  shown  in  Equation  (2)  below. 

FTCHzT  - 2h  /E^  % PT5HiT 


Figure  8 shows  that  this  is  a fairly  good  approximation. 
Dielectric  constant  scaling  results  from  an  approximation  in  the 
second  term  of  Equation  (1)  as  shown  in  Equation  (3). 


( 


1 + .61 


(E^-1) 


(B/H) 


.1155 


\l/2 
) * 


(B/H  « 72) 


(3) 


Figure  9 shows  that  this  is  also  a fairly  good  approximation. 

The  basic  design  tool  becomes  then,  a graph  of  Equation  (1)  with  the 
resonant  frequency  plotted  against  the  A/B  ratio,  with  the  pareuneters 
'A',  'H',  and  'Er'  specified  (graphs  where  A = 3.7,  Er  = 2.5,  and 
H .015,  .031,  .062,  and  .125  are  shown  in  Figure  10).  Once  this 
graph  is  used,  frequency  scaling  and  dielectric  constant  scaling  are 
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employed  to  arrive  at  a close  solution.  The  solution  for  'B'  can 
then  be  plugged  back  into  Equation  (1)  to  remove  any  errors  that 
might  have  resulted  from  the  approximations  and/or  from  reading  the 
graph. 

Design  Procedure  for  Size  Calculation 


Step  1.  Given 

a.  Substrate  thickness  = 'H* 

b.  Desired  resonant  frequency  E '^o' 

c.  Relative  dielectric  constant  of  the  substrate  e 'e  ' 

r 

Step  2 . Select  an  A/B  ratio 

Step  3.  Read  the  frequency  'F'  for  the  appropriate  substrate 
thickness  in  Figure  10 


Step  4.  Multiply  'F'  by 


Step  5.  Calculate  A = 3.7  (F/F^) 

Step  6.  Calculate  B = A/ (A/B  ratio) 

Step  7.  To  remove  any  errors,  substitute  'B'  into  Equation  (1)  and 
calculate  'A' 


FREOUENC'f  (GHZ) 
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Flu.  A-2  LENGTH  "A"  VERSUS  THE  SQUARE  ROOT  OF  Er  WHERE : 
B = A,  FREQUENCY  = 1 GHZ,  H « 0.047  INCHES 
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:SIGN  CURVES  FOR  VARIOUS  SUBSTRATE  HEIGHTS  SHOWING  THE  RESONANT  FREQUENCY  VERSUS  THE 
lO.  WITH  THE  PATCH  LENGTH  "A"  FIXED  AT  3.7  INCHES,  AND  THE  DIELECTRIC  CONSTANT  Er  FIXED  AT  2JB 
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Fine  Tuning  the  Calculation 

When  a particular  batch  of  material  is  going  to  be  used  to 
construct  several  antennas  of  differing  frequencies,  the  dielectric 
constant  used  in  Equation  (1)  can  be  adjusted  such  that  it  "predicts" 
the  exact  frequency  for  a test  antenna.  This  is  done  by  solving  for 
Ej.  in  Equation  (1)  (shown  below)  and  then  using  the  parameters  and 
the  measured  frequency  of  the  test  antenna  to  calculate  an  effective 
dielectric  constant  for  the  substrate. 


W = .61* 

V = (F*H*10"^)  ^ 

U = (A*F*10“^)^ 

X = 34.82684258  + (W-1)*V 

Z = (4*V  - U*W) 

Y - X*Z  -278.6147406*V 

2 r 2 2 2 2I 

E^  = -Y/Z^  - (Y/Z^)  ^ - (X/Z^)  ^1 

We  are  investigating  the  accuracy  of  this  procedure  in 
determining  the  dielectric  constant  of  a substrate. 

Input  Impedance  Calculation 

The  calculation  of  the  input  impedance  at  resonance  given  in 
reference  (3)  proved  to  be  inaccurate.  However,  knowing  the 
impedance  at  a given  feed  point,  the  impedance  at  any  other  feed 
point  can  be  determined  using  Equation  (4) . 


R.  = R sin  (Y  *P) 
in  o ' o 


(4) 


where  P = 

A + 2H/  /e~ 


and  Ro  is  found  experimentally. 
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The  design  procedure  would  thus  be  as  follows: 

Step  1.  Determine  desired  input  resistance  based  on  system 
considerations 

Step  2.  Build  a microstrip  antenna  at  the  appropriate 
frequency  and  feed  it  at  an  arbitrary 

Step  3.  Measure  the  input  resistance 

2 

Step  4.  From  Equation  (4)  calculate  = R^/sin  (Yj^*P) 

Step  5.  Calculate  'Y  ' for  the  desired  input  resistance 
using  EquatiSn  (5). 
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Function 
All  other 


APPENDIX  B 

Computer  Programs  and  Subroutines 


Programs 

1.  A OF  B 

2.  A OF  ER 

3.  A VS  F 

4.  DRIVE 

5.  E AND  H 

6.  F VS  ATOB 

7.  MICROAN 

8.  PATTERN 

9.  RESULTS 

10.  RIN  VS  H 


Subroutines 

1.  ANGL 

2.  CONSTAN 

3.  DIELECT 

4.  LENGTH 

5.  PAT 

6.  RIN 

7.  SIMP 
8 , T 

9.  U 

10.  WATTS 

11.  XFIL 


XFIL"  was  taken  from  the  NSWC  user  library  and  modified, 
programs  were  written  by  John  McCorkle. 
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PROGRAM  AOFe(OUTPUT, INPUT) 

C THIS  program  plots  LENGTH  A VERSUS  WIDTH  B 
DIMENSION  EC(A),HC(.T),AL(51),RX(51) 

DATA  EC/?.,?,?5')n07e-29,?.55,?.‘314?4  730  3/ 

DATA  HC/.  03175,  ,046B7F,  . (1625/ 

A (H  ,ER,F,9)  = ( (1 ,18F10-P*4,»H»SQFT (ER) )/ (2 , * F* SQRT ( 1 . ♦ . 6 1* (ER-1 .) ♦ 

1 (“/H)**. 1155)) » 

F = 1 .E9 
FR=FC(3) 

DO  12  J=l,3 

20  FORMAT(/,»  A P*) 

PRINT  20 
H=HC(  J) 

B = 1 . 

no  10  1=1,51 
AL( I) =A (H,ER,F, P) 

9X( I) =P 

PRINT  2,AL(I)  ,RX(I) 

2 FORMAK  2 (lPG2n.  5)  ) 

10  R=P+.15 

IF(J.LT.2)  CALL  C ALC M 1 ( 51 , B X , AL , 0 , 1 . 0 ,8 .5 , 3 . Ha , 3 . 75 , 12, 5 , 9. 0 , 
119HMICPCSTRIP  ANTENNAS, -19, 15HL5NGTH  B INCHES, 15, 

2 15HLENGTH  A INCHES,  1 5 , 1,  18  ) 

IF<J.EQ.2>  CALL  CALCM 1 ( -c j , By, AL , -1 ) 

IF(J,E0.3)  CALL  f ALC H l ( -5 1 , BX, AL , -2 ) 

IF(J.GT.3)  CALL  C ALC M i ( -? i , BX , AL , -5 ) 

12  CONTINUE 

CALL  GPID(0.,0. ,.5,.5,25,  18) 

CALL  CALCM1(0,0.) 

PND 
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PROGRAM  AOFER  (OUTPUT ,TAPEe=OUTPUT ) 

EXTFRNAL  PUN 
COMMON  H,ER,F 

OIMENSICN  CONET  C^)  , AL  (51)  ,RA  (51  > , ATOBd  > , A ( 14>  , ANS  (1  > 

NAMELIST  /PARAM/H,r, ATOP 

DATA  CONST/. 031?E,. 046875 0625/ 

OAT  A A/1 .,6.,1 . E-1,1  . l.E-4, .5 5»1.  Ee,q9.0»0. ,1 .,1.  1, 20./ 

F=i  .pq 
EINC= .144 
H=T./64. 

J = 1 

ATOB( 1) =1  . 

EP=?. 

WRTTP  (6,1» 

1 FOPMAT(/) 

WRITE (6 ,PARAM) 

WRITE  (6,2) 

PORMAT(  5X,*FP*,10X,*A*,7X,*1/SOPT  (£»)♦) 

POPMAK  3(510.7)  ) 

00  10  1=1,51 
BA (T)=1./SGRT(EP) 

CALL  ROOTER  ( FUN, A, ATCO  ,AMS ) 

AL(I)=ANS(l)  •ATOB(I) 

WRITE (6 ,3)ER, AL  (I) ,BA  (I) 

10  ER=ER+EINC 
Xl= .7071068 
X2= .3296902 

IF(J.LT.7)  CALL  CALCM1(51,RA,AL,0,X1  ,X2  , 1 . 5 , 4 . 2 , 1 2. 0 , 9 . 0 , 1 9HM ICR 
lOSTRIP  ANTENNAS, -19, 1 OMl/SQPT (£R) , -9, 15HLENGTH  A INCHE S ,1 5 , 1 , 1 8 ) 
TF(J.E0.2)  call  CALCM1(-p1,PA,AL,-1) 

IF(J.EQ.3)  CALL  C ALC^  1 ( -«;  1 , P A , AL  , -2  ) 

TF(J.GT.3)  call  CALCH1(-51,PA,AL,-5) 

12  CONTINUE 

CALL  GRin(0.,0. ,.5,.5,24, 18) 

CALL  CALCM1(0,0.) 

ENO 

FUNCTION  FUN(P,CUMY) 

COMMON  H,PR,F 
dimension  OUMY(l) 

9N=B*0UMV(1) 

FUN=( (1 .16E10-F*4.»H*SQRT (EP) )/ (2.*F*SQRT (1  .♦.SI* (ER-1 . )♦ (B/H)** 
1.1155)))  -BN 
RETURN 
ENO 


NSWC/WOL/TR  76-69 


PROG*»AM  AVSF  (OL/fpuT,  INPUT) 

C***THIS  program  PLOT'  AND  PRINTS  <A<  VFRSES  FpEOUcNCY  AT  A SPECIFIED  A/B  PATIO 
c***ano  various  values  of  H. 

C***INPUT  DATA  IS  - A TOB,F M IN ,FMA X , NP 

C***NP=NUMRER  OF  FPEQLENCI'S  (POINTS)  TO  BE  CALCULATED. 

C***  NOTE  (F)'AX-FMIN)/{NP-1(  = CONVENIENT  INCREMENT  FOR  THE  FREQUENCY 
C 

COMHON/e/A,P,F, H, EP 

DIMENSICN  HH(<4)tAl(?n,A?t51),A3(?l),A4(51),FX(51) 

DATA  HH/,01E,.OTl,.Dfi?,.l?5/ 

DATA  Xl,X2,Yl,Y?,XL,YL/,l,10.,t.,10.,-26.,-9./ 

M0D(I)  = I-(1/4)*i. 

HH< 1) =. 047 
ER=2.  5 

READ*  ,ATOS,FMIN,FMAX,NO 
FINC= (FMAX-FMIN)/  (NP-  1) 

B = n . 

FsPMIN 

PRINT  1 ,ATOP,EP .HH(1 ) ,HH( 2)  ,HH(3) ,HH(4) 

1 FORMAT(  2X,*AT0R  = * ,F6. 3,  EX  , • EF  = * , F 6 . 3 , / , 22X , 4 ( 1 X , ♦Hr*  , F5 . 3)  ,/, 

I 14X,*FRED(GHZ) •,4(?4X,*A*) ,/) 

DO  10  1=1, NF 
FX( I) =F*1 .E-3 
Hr  HH  (1  ) 

CALL  LENGTH(ATnP,?0) 

Al( I) =A 
H=  HH (?) 

CALL  LENGTH (ATOP, 20) 

A?( I) =A 
H=  HH(3) 

CALL  LENGTH(ATOT,?0) 

A3( I)=A 
H=  HH ( 4 ) 

call  LFNGTH(ATn?,20) 

A4( I) =A 

PRINT  ?,FX(I),A1(I),A?(I),A3(I),A4(II 

2 FORMAT(2X,5 (1PG3S.8> ) 

IF(M00(  I)  .EO,  (I)  PRINT  14?? 

1425  FQPMATdH  ) 

10  F=F*F1NC 

call  CALCMl (NP,FX,A1 ,0, XI ,X2,Y1,Y2, XL,YL, 

119HMICRCSTRIP  antennas, -19, ISHFREOUENCY  (GHZ), -15 
?,15HLENGTH  A IN CHfS , t E , 1 , 18 ) 

CALL  LAeLQG(?8. ,.1,1Q.,0.» 

CALL  LABLOG(9.,l.,10.,D0.) 

CALL  GRID(0.,0. ,14. ,9. ,-?,-!) 
call  CALCM1(0,0.) 

END 
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C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


4 

1 


5 


15 


?0 


PROGPAM  ORIV^' (INPUT, OLTFUT) 

THIS  PROGRAM  PRODUCES  PLOTS  OF  R-I NPUT , EFF IC IE NCY , AND 
PANOHIOTh  versus  the  parameter  chosen  by  tOP^t 
OPT=T  FOP  PLOTS  . VS . A/P  RATIO  AND  INPUT  DATA  IS. 

F,  H,  ER,  XYO,  LOSSTAN,  ABMIN,  ABMA X,  NPOINTS. 

OPT=1  FOP  PLOTS  VS.  FRECUENCY  ANC  INPUT  DATA  IS 
A/P,  XYO,  H,  -K,  LOSSTAN,  FMIN,  FMAX,  NPCINTS 
OPT  = ‘’  FOR  PLOTS  .VS.  ER  AND  INPUT  DATA  IS. 

F,  H,  A/B,  XYO,  LOSSTAN,  ERHIN,  ERMAX,  NPOINTS. 

OPT=^  FOR  PLOTS  .VS.  RR  AND  INPUT  DATA  IS. 

A,  P,  YO,  LOSSTAN,  RRHIN,  ERMAX,  NPCINTS. 

REAL  LAMOA,LAHOAr,,L,LCSSTAN,IM 
COMMON/P/A ,P,F, H,ER 

COMMON/E/RS,RC, P I ?, I H ,L , PI , L AMD  A G , L AHOA , LOSS! AN , ZO , PA , PP , P 

common/f/eff, oeltaf, gain 

DIMENSION  X(15n),Y(150),D(15n),E(150) 

DATA  FREO/PHFPEOUENCY/, APR/9HA/e  P A T 10/ ,FeR/gHE PSI L ON  R/ 

OAT  A IM,A0INC,FINC,EPINC, Y1 ,Y?,XL,YL/1.,0.,0.,0.,0.,0.,12.0,9.C/ 
DATA  TANG/PHLOSS  TANG/ 

FREQU  (H,ER,  A,P)  =5.gOt4‘»7165£9/(A*SQRT(l.*.6l*(ER-l)*(B/H)**.ll55) 
1+2. •H*SORT(FR)) 

MOD (I)=I-(I/4)*4 
CALL  ANGL 
READ*  ,NT 
PEAD* ,OPT 

if(opt.eo.-ioo)  go  TC  4 
A0TNC  = O . 

FINC=0. 

FRINC  = 0 . 

IF(OPT.EO.l)  GO  TO  5 
IF( OPT. E0.2.)  GC  TO  IE 
IF(0PT.E0,3.)G0  TO  20 
IF (OPT.LT.O. » GO  in  12C 
REA  0» ,F ,H, E P, XY C, LOSS  I A N, AP MIN, A0MAX, NPOINTS 
NAMELIST/CPTO/F ,H,ER, LOSS  tan, XYC  ,NT 
AX=  APR 
Xl=  AOMIN 
X2=  ABMAX 

ABINC=(X?-X1I  /NFOIMS 
AnB=Xl 

P = 5 ,75Eg/ (F*ATOP*SOPT  (FR) » 

PRINT  OPTO 
GO  TO  400 

READ* , A TOfl, XYO, H,ER, LCS STAN, FHIN,FM AX, NPOINTS 

AX=FRErj 

X1  = F'«1N 

X?=FMAX 

FINC- (X2-X1)/NP0INTS 
F = X1 

P=5.75Fg/(F*AT0o*S0RT(RR) ) 

NAMELIST/OPT1/ATO0,XYC,H,ER,LOSSTAN,NT 
PRINT  OPTl 
GO  TO  400 

READ* ,F,H,ATOe, XYO, LOSSTA N, EPMI N , ERMAX , NPOI NTS 

AX=EER 

X1=EPMIN 

X2=fpmAX 

FR=X1 

0=5 .7  5Eg/ (F* ATOP* SORT  (fr) ) 

ERINC=(X?-X1) /NPOINTS 

NAMELIST/CPT2/F ,H , AT  OP , XYO , LOSST A N, NT 
PRINT  OPT? 

GO  TO  400 

READ* , A, B,H,YO, LOSST AN, EPMI N, ERMAX, NPOI NTS 
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flTOP=  A/P 
AX=E- R 

xi=erhin 

X?=ERHAX 

ER=X1 

ERINC=( X2- XI » /N POINTS 
NAHELIST/OPTT/A ,B ,H, LGSSTAN , Y0,KT 
PRINT  OPTT 
400  N=NPOINTS+l 

456  PO=MAT(/,10X,»A*,14X,*A/?*,15X,»F*,14X,*ER*,  9X, ♦DELTA  F* , 6X, 

l*EFFICIENCY»,tlX,*  70  ♦,  8X,*R  INPUT*,/) 

PRINT  U56 
DO  100  1=1, K 

IF( 0PT,EQ.3.)  FiFRFOU (H,FP, A,B) 

TF(0PT.EQ.3.)  GO  TO  10 
GALL  LENGTH(ATOe,20) 

YO=  A/ XYO 

10  CALL  CONSTAK(3. 19E-8 , 1.47E6) 

CALL  WATTS(NT,M) 

CALL  RIN(-1.,W, RH) 

F (I )=EFF 
0(1 )=OELTAF 
FG=P*  1.  E-<3 

FIX=( 4. 1316e667E-3/Hfl.5822-H*7.466026667)/ ( FG** ( . 3281 4+2 . 7 515* H+ 
1 6.  90  77*H*H)) 

RI=RH*SIN(YC*P)  **2  *FIX 
Y (T  ) = RI 
X (I  ) = ER 

IF( OPT.EQ.l.)  X(I)=F 
TF( OPT. £0.0  .)  X(I)=AT09 
789  F0PMAT(8(lPGie.  5)  ) 

PRINT  7 89,A,AT0P,F,ER,0(T),E(I)  ,ZO,Y(I) 

TF(MOO(I)  .EC.C)  PRINT  ? 

2 F0'?NAT(1H  ) 

F=e+PINC 
ER=  ER  + ERINC 
100  ATOB= ATOB+APINC 

CALL  CALCHKN  , X ,Y , - 1 , XI , X2 , VI  , Y2  , XL  , Y L , 1 BHI-ICROST  RIP  ANTENNAS, 

1 -19, AX  ,-9,3HRIN,3,1,18) 

CALL  GRin(0.,C.,.5,.5,?4, 19) 

CALL  CALCMKOjO.) 

GO  TO  1 

CALL  CALCHKN  , X,0 , - 1 , XI , X2,  Y1 , Y2  , XL  ,Y  L , 19HH ICPOST  RIP  ANTENNAS, 

1 -19,AX  ,-9,7H0£LTA  F, 7, 1,18) 

CALL  GRIO(0., 0. ,.5,. 5,24, 18) 

CALL  CALCHKN  ,X  ,E  , - 1 , XI , X2,  Y1 , Y2,  XL  ,Y  L , 19HHICR0ST  RIP  ANTENNAS, 

1 -19, AX  ,-9,inHEFFiriENCY,10,l, 18) 

CALL  GRIO(0.,(J.  ,.5,.5,24,  18) 

120  STOP 
fNO 
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OfclOCPAH  - ANOH  (INPUTtnjTPUH 

C*»»H»SC0  ON  NO?  K4L0I<S  PA^tR  ON  PtCROSTRIP  ANTENNAS  1975 

C»*»ThIS  PROr.RAH  SNOWS  HOW  CURCUIT  BOARO  TOERANCES.  VARIATIONS  IN  ER  ANO  H, 
C»»*WIH.  fFFtCT  TNE  operation  OF  THE  MICROSTRIP  PATCH  OUTPUT  DATA  IS 
r •••FP  eg  (C|  -.WEIGHT  <H) --DIELECTRIC  CONSTANT  (tR»  — INPUT  RESISTANCE  |RI»  —GAIN  — 
r»»»FFFICIFNCT(EFFI  — PA NOHIO TH (DEL T AF » --TO 
r»**INPUT  OATA  IS  AS  FOLLOWS. 

C»**4 ,e,yO,L0SST4N,MHlN.MNAX,NH,ERNlN.ERNAX,NE.KK 

r*.»TMr  ppogfan  continues  reading  new  input  cards  until  none  are  left 

r •••NH-N'IWPER  pF  VALUFS  (POINTS)  FOR  H 
r»»»NE-Nt)MnrR  OF  VALUES  (POINTS)  FOR  ER  TO  TAKE  ON. 
r»»»  AOTf.  ( M«AX-HMIN) /( NH- 1)  =CONVIENIENT  INCREMENT  FOR  H 
(;•••  (EPMAX-:pMIN)/(N'--1)  = rONVENIENT  INCREMENT  FOR  ER 

r»»*KK=i  Fr?  A SERIFS  OF  H valufs  retheen  : ach  Change  in  er 

r»**KKe?  FOR  A SERIFS  OF  "R  VALUfS  BETWEEN  EACH  CHANGE  IN  H. 

C 

RFAL  K. LAMOA.LAMDAG.mJ.L.LOSSTAN, IM 
COhmon/o/a.b.e.h.ep 

CONMON/ '/RS ,RC, PI2, IH.L.PI.LAMOAG.LAROA.LOSSTAN  ,70,PA,PP,P 
CONMON/F/FFF.nF  LTAF .GAIN 

FRECIKH.fr,  A,9>  rS.O01i.271b5E9/)  A*SGRT  (cp.i)  . ,11551 

1*?.»h»S0°T(FR)» 

IH=1. 

CALL  ANSI 

1 OEA0*.A.0,y0.L0SSTAN.HMIN,HMAX.NW  , ? R M IN, E R MA X. NE  . KK 

IF<roF(SLINPUT) .NE.C)  GO  TO  °9 
NAN, L IS  T/PAPAM/ A, 0, TO. LOSSTAN 
rR=FPNiN 
H=HMIN 

riNC=  (er'^ax-ernin)/  (NE-i) 

MINT:  (HNAX-HMIfU/INM-l) 

PRINT  piPiw 

PRINT 

IF(KK.E0.2)  go  TO  1,0 
DO  i:  ItI.NF 
h=hhin 

00  ?C  Jel.NH 
F = EP20ll(H.tR.A.P) 
call  CONSTANO.  lOE-A.l  .RTEb) 
call  WATTSd.W) 

CALL  RIN(-1 ..W.PH) 

PRINT  1 .F.H.ER.RH.GAIN.fFF.OZLTAF.VR.YF 
?(!  H=H*HINC 

PRINT  2 

IP  ER=fr,£TnC 

GO  TO  I 

1.0  DO  5C  1 = 1 .NH 

ER=FRMIN 
PO  EC  J=l.Nc 
F=FPFQll  (H.EP,  A.  P) 
call  CONSTANO.  19E-A,1  .I.7E6) 
call  WATTS(I.H) 

CALL  RIN(-1 ..W.RH) 

TE  = RH»( ST N( (Y0».01) ‘P) ••2-SIN ( ( Y0-. C 1 ) •P) **  2 ) 

VR* AS IN (SORT (50  ./RH))/P 

PRINT  3 .f.H.fr, PH.GAIN.FFF.nFLTAF.VR.TF 

EC  FR=FP*fINC 

PRINT  2 

50  H=H»HINC 

GO  TO  1 
FORMATdH  ) 

FORMATOdPGlW.F)) 

FOPMAT( 3X,»FPEQUENCY*,  9X ,• H* , 1 IX ,• ER • , 11 X , • R *,  8X,*GAIN^. 

IflX, •EFFICIENCY*. 5X,^BANOWIDTH*,3X, •SO  OHM  YO* , 3X,*. 01 XX  ERROR*,//) 


99  STOP 
END 
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jgOM  OOP!  fJJBtilSiffiiJ  TP  


PROGRAM  FVSATOn  ( OUTPUT NPUT ) 

C***THIS  PROGRAH  PLOT?  ANO  PRINTS  THE  RESONANT  FREQUINCY  VERSES  A/0  RATIO 
C***F0R  a fixed  value  of  a and  VARIOUS  VALUES  OF  H. 

C***INPUT  DATA  IS  AS  FOLLOWS. 

C***A,ER, A0MIN,AflHAX,NP,Yl, YE 

C***NP=NUM0ER  OF  A/P  RATIOS  (POINTS)  TO  OE  CALCULATED  AND  PLOTTED 

C***  NOTE.  ( ABMAY-ABMIN)/ (KP-1)  = CONVENIENT  INCREMENT  FOR  THE  A/0  RATIO  POINTS 

C 

dimen  SIGN  X(51)  ,F1(E1),F?(51),F3(51),F4(51) 

DATA  HI ,H?,H3,HA/.ni5 ,. 031, .062, . 125/ 

MOO  ri)= I-(I/5)*K 

FREQU(H,ER,A,P) =5.90 1a27165E9/(A*SOPT(1 .♦.61*(:R-1)*(B/H)**.1155)* 

1 2. •H*SORT(FP)) 

100  READ* ,A,rR,AeMIN,ABMAX,NP,Yl,y2 
IF(E0F(5LINFUT>  .NE. 0)  GC  TO  999 
PRINT  1,A,EP,H1,H2,H3,HA 

1 FORMA T ( 10 X,*A  = * , F6. 3 , 5 Y ,* ER  = * , F6 . 3 , // , 1 f X ,♦ AT  09* , 3X , 4 ( 1 EX 
1 ,*F=>En  (GHZ)»),/,21X,4(15X,*H  = *,F5.3)  ,/) 

AINC= (A9MAX-ABMIN)/(NF-1) 

ATOB=  ABMIN 
DO  10  1=1, NP 
RrA/ATOP 
X(T )= ATOP 

Fl( I) =FPEOU  (HI, CD,4, n) *1. t-g 
F2( I) =FREOU (H2, -P,A, P) • 1. E-9 
F3( I) =FPEQU(H3, EP,A,n)*l.E-9 
F4( I) =FREOU (H4, EP,A, R| * 1. E-9 
PRINT  2,X(I),F1(I),p2(I),F3(I),F4(I) 

2 FOpMAT(2X,5 (1PG25,9) ) 

IF(MOD( I)  .EO.O)  PRINT  14^5 

1425  FORMATdH  ) 

10  AT08= ATCP+AINC 

CALL  CALCMl  (NP,  X,F1,  0 ,APMIn,A8m.''X,V1,Y?  ,12.  ,9., 

1 19HMICR0STRIP  ANT-nnAS  ,-19,9HA/e  R AT  10  , 9 , 9HFR EQ'JENCY , 9 , 1 , 1 8) 

CALL  CALCMl  (-NP,X,F7,-1  ) 

CALL  CALCMl (-NP,X,FX, -2) 

CALL  CALCMl (-NP,X,F4, -5) 

CALL  GPin(0.  ,0 . , .F, i 5, 24,1 8) 

CALL  CALCMl(0,n.) 

GO  TO  100 
999  STOP 
PND 
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program  ••rCROAN  lINPlJT,  OUTpJti 
C TmIS  PR0GRA*1  designs  « BACK  F$0  HXCROSTRIP  ANTENNA. 

C INPUT  DATA  IS  AS  FOLLOWS, 
f DESIRFD  FSrOUENCV  IN  NM7 
C A/e  RATIO 

C SU«STRATE  M«IGNT  IN  INCMiS 
C OTELECTPIC  CONSTANT 
C OESIRFO  input  impedance 

C HEASUREP/CALCULATEO  INPUT  RESISTANCE  FOR  A TEST  ANTENNA  <MOPEFULLT  l.OT 
C LOSS  TANGENT. 

r THE  PROGRAM  IS  S-T  UP  TO  RUN  ON  POTm  MiCH  SPEED  LINE  PRINTERS  AS  WELL  AS 
C TFLETyPf,  FOR  LINF  PRINTER  USE.  TMf  FIRST  FREOUENCT  READ  IN  SHOULD  BE 
C N-GATIVF.  TMF  PROGRAM  CONTINUES  READING  UNTIL  IT  RUNS  OUT  OF  DATA  CAROS. 

C 

REAL  K. LAMOA.LAMOAG.HJ.L.LOSSTAN. IN 
COMMON/ P/ A, n.F.M.ER 

C0HMON/=-/Rs,cc.pi2.  im.l.pi.lamdac.lamoa.losstan.zo.pa.pp.p 

COMMON/F /erF, DEL TAF .GAIN 

M0n<llsl-ll/4l*<. 

iMsl, 

call  ANGL 

NsO 

PRINT  eo 

TO 

1 REAP* ,C .ATOP.H. fR.ROS. * .LOSS! AN 

IP  (C  OF  (Ct  INPUT!  .NE.-J*  GO  TO  !»R 
IsT*l 
FG  = E 

TF(F,GT0,I  go  to  I-J 
Ns  I 

PRINT  «»: 

FGs»Pr, 

I*;  FsFf,*l.PT 

CALL  L.Nr.TwfATOR.E')! 

TOsA/ 3, 

MsM*. O'  I 

call  CONSTAN{1,l<)E-^.l*‘»Trf‘» 

CALL  WArrSU.M) 
call  RIN(-1 

TIT*  tW.  I 31f*f6hrP*  3/H*l  CFG*  •!.  T78lA»2.r51S*M* 

I 

PlsP 

PMIsOm*  FI t 

MSM-. 3C  I 

call  const an  I. 3. IRE-A.l , 47Ff » 

call  mattsu.wi 

call  RI Nl-l , , M, Rm» 

FITS  14.  I - 3/H*l  .^^?2-m*7.‘.«.6C?6667»/  CTG»  •<.3?fll4.*2.75l5*M^ 

1 f»,  RC  »7"m*m>  I 
RMspH»PTv  •le 
rOsA/?, 

7POsSORTcROS/Rm> 

IF  ( 7PO.G’^.l  .1  RnSsRM*<  STN<TO*Pl  ••  ?l 
IF  f 7*>0.  L’.l  . > vns4SI»  I 7PO»/P 
TlsA/2,-vO 

YO“  Fs  R:)^-RM»<:  IN  I I TO-.  3 1 » ‘P*  ? 

M-osCMj  •<Pi^<(ro*Pt»  -RCS 
pPLTAFsn  L^AF*!  ,r-9 

IF<M.Lr,l|OFINT  A.T.TO.Tl.VOe-.HFt'.O'^LTAF.r.a  tN.RnS.FG.ATOe.M 

I .'■o.LOSS^AM,*  ,RM  ,o 

IF  IN.r.T  .-j » Pi  IN  T 30,4,  R,m«R  ,*<,cR,Fr.,t  OS,V  1,  Rh,VO:R,LOSST  AN,  X, 

1 0-lTAF,  GAIN.EFF 
IFfN.LT.ll  PCINT  72 
IF CMOni I» ,?c. 01  PRINT  SO 


GO  TO  1 

3C  FORMATCl« (1I,F8,4»» 

LC  FOPMATf»l  A*,  SK,*M  ' PRO*^  • • GX  , 'H*  , F X , • E »•  . S T , •F  • , SX  , • R OS*  , 

lT*,*Yl*,*X,»B0*,AX,*yO  I PROP*  ,3X**L0ESTAN*, St,*  t*  , JX,*eJ NOWlOTM  • , 
l*f.ATNIO‘3l*.3X,*'’FF*| 

«C  FORMAT!*  *1 

PC  POPRAT!*  TVPf  IN  Tm£  FOLLOWING  PA  RA  •»'  T^RS , * , / , • FRED  IN  GM7*,/, 

1*  A/n  RATIO*,/, • HtlOMT  IN  INCmE<^*./.*  OIELFCTRIC  CONSTANT*,/, 

?•  nrsiRPO  INPUT  OESIST  ANCe*  */••  '"ORR'CTION  FACTOR  * “FAfOREO/*. 
l*rAirULATFO  INPUT  RESISTANCE  OF  A T^'ST  ANTENNA*  ,/, 
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1/.*  Y0».T3Q,***,Fl2.f,/,*  Yl»A/2  • Y" • , T 3C . . Fl 2, f ,/, 

?•  OMMS  P.'R  .'*15  CHANGE  IN  YG  • , T J 1 , • s • , F I , 

3*  OMMS  »“R  .OOITX  CHANGE  IN  M*  , T3  3 , • » * .Fl  7,  f>  ,/, 

4*  -300  TANOwTOTm  IGM7» * .T30 ,F12, -,/,*  GAIN  !nO  OVER  ISOTOOPtCl 
R».T30.*8*,Fi?.f,,/,»  input  RTRlSTANC'*  .T30,*s*,Fi7.f.,/, 

A*  FRFO'J-nCY  IN  r,M7*  ,T30  .*«*  .F  12,F  ./ »•  A/0  P A T 10  * . T * 0 . , F I 2, 6. /, 

T*  SUPSTSATP  mpiGmT  •.T  If  ,*s*,Fl  7,‘>, /,*  OIELtCTPIC  CONSTANT* 
O.TSC,*i*.fl2.f>,E.*  LOSS  TANGf  NT*  .Tjr  ,•*•  ,f  I 2,F,, /,•  fORRfCTlON* 
0*FArT0R*,T3C.*»*.Fl2,e./.*  RC  * • T3a , • • .Fi  ? , t • /,  * P*  , T3C  • .^12  .f>  I 

OR  STOP 
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Prograa  to  Design  Microstrip  Antennss 
With  a Back  Feed 


By:  John  W.  McCorkle. 

IS/It/IP.  I0.t5.2a* 

IWIIE  3AK  A5PC/397 
USEFNUMFF.h:  aab. 

lEKMINALt  53. T1Y 

hECOV/Eh  /SYSTEMtbAICM 
ShFL.POCOC. 


hFL. aOOCO. 

/'C-E1.MlCh.'flN/L.M.6q5. 

/CET.Lli-^y>P:r/L'J=^=<5. 

/LINK. F=.v!Uf<3AM.H. LIE. X. 

TYEE  IN  7KE  fOLLJEINC  EAFAMETEAS. 

FFEC  IN  CHZ 

A/t  FATIO 

HEICm  IN  INCHES 

DIELECTMC  CONSTANT 

DESIFEO  INEDT  FESTANCF 

COhFLCTION  FACTOh  (USUALLY  I.O) 

LOSS  TANCENT 

T 1.2  1.25  .125  2.5  SO.  I.  .0005 


length  a 

m 

2.990736 

WIDTH  b 

m 

2.392589 

YO 

m 

.383353 

Yl*A/2  - YO 

m 

I . 1 1 201  a 

OHMS  FEF  -Ol"  CHANCE 

IN  YO  « 

2. aT99aO 

OHMS  FEh  .cot”  CHANGE 

IM  M > 

-.295700 

-3Db  bANUH  DTH  < GHZ) 

9 

.028777 

CAIN  (IE) 

m 

5.29a9a7 

INFUT  ffsi stance 

m 

50.000000 

FFFOUENCY  (GHZ) 

m 

1.200000 

A/b  FATIO 

• 

1.250000 

SUbSTFATF  HEIGHT  H 

m 

. 125000 

CIILECTMC  CONSTANT 

m 

2. 500000 

LOSS  TANGFNT 

■ 

.000500 

COKFECTION  FACT3F 

■ 

1.000000 

AH 

■ 

35A.968a.63 

P 

■ 

.997695 

IF  YOU  HAVE  N3  HOFF  DATA  TYFF 

IN  'STOP* 

7 3.  1.1  -OAT  2.5  50.  1 

. .001 

LENGTH  A 

s 

1 . 1951  78 

WIDTH  B 

• 

1 .086526 

YO 

m 

. 1 a9  638 

YI«A/2  - YO 

m 

. aa795l 

OHMS  PEF  .Of  CHANGE 

IN  YO  ■ 

6. 36aa03 

OHMS  PEh  .001 ••  CHANCE 

IN  H » 

.336172 

-3Cb  EANDV.IDTH  (GHZ) 

w 

.077330 

GAIN  (bfc) 

m 

5.  150967 

INPUT  FFM  STANCE 

m 

50.000000 

FPECl'ENCY  (GHZ) 

9 

3.C00CCG 

A/P  KATI3 

m 

I.IOOCCO 

SUESTFATF  HEIGHT  H 

m 

.OaTOCO 

DIELELTKIC  CONSTANT 

■ 

2.5GCCC0 

LOSS  TANG!  NT 

• 

.001000 

COFFECTION  EACTOF 

• 

1 .OfCCOO 

FH 

• 

373.280678 

P 

M 

2. 506001 

IF  YOU  HAVF  NO  M.tFF  DATA  TYPE 

IN  'STOF* 

» AUt: 

•lEFMlNATFD* 

AEFFJX  COM  OF 
FETIFNI TZi/CJ) 


nv  IS  1 


. 57 


M5 

M5 


LUC  OFF.  10. 31 • /A. 
CF-  3.flHA  EEC. 
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X-Y  PLANE 


PHI' 

EPHI 

ETHETA 

ETOTAL 

scale  40  00  = 

5.06136 

-IflO.OO 

17.729 

-411.04 

17.729 

-175.00 

17.451 

-212.42 

17.451 

-170.00 

16.562 

-206.45 

16.562 

-165.00 

14.867 

-203.01 

14.867 

-160.00 

11.816 

-200.63 

11.816 

-155.00 

5.2929 

-198.85 

5.2929 

-150.00 

-4.4025 

-197.44 

-4.4025 

-145. CO 

10.601 

-196.32 

10.601 

-140. CO 

16.231 

-195.40 

16.231 

-135.00 

19.854 

-194.64 

19.854 

-130.00 

22.548 

-194.02 

22.548 

-125.00 

24.695 

-193.50 

24.695 

-120. OC 

26.476 

-193.09 

26.476 

-115.00 

27.994 

-192.75 

27.994 

-110.00 

29.313 

-192 .49 

29.313 

-105.00 

30.474 

-192.29 

30.474 

-100. GO 

31.509 

-192.15 

31.509 

-95.000 

32.440 

-192.07 

32.440 

-90.000 

33.283 

-192.05 

33.283 

-05.000 

34.052 

-192.07 

34.052 

-80.000 

34.756 

-192.15 

34.756 

-75.000 

35.403 

-192.29 

35.403 

-70.0CC 

35.998 

-192.49 

35.908 

-65.000 

36.547 

-192.75 

36.547 

-60.000 

37.052 

-193.09 

37.092 

-55.000 

37.516 

-193.50 

37.516 

-50.000 

37.941 

-194.02 

37.941 

-45.00C 

38.326 

-194.64 

38.326 

-40.000 

38.672 

-195.40 

38.672 

-35.000 

38.979 

-196.32 

38.979 

-30.000 

39.247 

-197.44 

39.247 

-25.000 

39.475 

-198.85 

39.476 

-20.000 

39.663 

-200.63 

39.663 

-15.000 

39.810 

-203.01 

39.810 

-10.000 

39.915 

-206.45 

39.915 

-5.0000 

39.979 

-212.42 

39.979 

0 . 

40.000 

0. 

40.000 

5.0000 

39.979 

-212.42 

39.979 

10.000 

39.915 

-206.45 

39.915 

15.000 

39.810 

-203.01 

39.810 

20.000 

39.663 

-200.63 

39.663 

25.000 

39.475 

-198.85 

39.475 

30.000 

39.247 

-197.44 

39.247 

35.000 

38.979 

-196.32 

38.979 

40.000 

38.672 

-195.40 

38.672 

45.000 

38.326 

-194.64 

38.326 

50.000 

37.941 

-194.02 

37.941 

55.000 

37.516 

-193.50 

37.516 

60.000 

37.052 

-193.09 

37.052 

65.000 

36.547 

-192.75 

36.547 

70.000 

35.998 

-192.49 

35.998 

75.000 

35.403 

-192.29 

35.403 

80.000 

34.756 

-192.15 

34.756 

85. one 

34.052 

-192.07 

34.052 

90.000 

33.283 

-192.05 

33.283 

95.000 

32.440 

-192.07 

32.440 

100.00 

31.509 

-192.15 

31.509 

105.00 

30.474 

-192.29 

30  .474 

110.00 

29.313 

-192.49 

29.313 
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PESULTSawWjTtOlTPUTl 

C THIS  PPOCPAH  COMPtLCS  *U  OF  THC  HE«SU*CO  DATA  INTO  A FORIt  MNICM  I.EMOS 
C ITSELF  TO  ANAITSTS. 

c input  data  is  as  fouous* 

C A,  «•  H , ~P»  HEASUPEO  FAEO  IN  NMl*  NEASUtCO  V0|  MEASUftEO  INHJT  RESTAMCE 
C ANO  an  iNTfGEP  JJ  HHEPF  CPOUPS  OF  JJ  LINES  ARE  FRIMTEO  MITN  A RLAMK  LINE 
C P<-TMFEN  EACH  GROUP. 

REAL  K(LAHOA,LANOAGtHL«L«LOSSTAN,lN 
COHHON/e/A«e»FtHtCR 

CONNON/c/eStRC»PIEitH*L»Pl»LANOAC»LAHOA«LOSSTAN»EO*FA»FR»F 
01HI--NSI0N  NFGn»»FORI19> 

HATA  p<fC/?h3m,?hpl»FHCU/ 

PRcrM'  |H,ER,A,Bt  sK,S01A?716SE9/|A«SORT(1.«.61*(£R>1>R|B/NI  R*.llf9l 

!♦?. •H»SOPT|F»l» 

Hoon,jMi-a/ji*j 

THSI. 

LOS?T*n«. 181 

PEAO  11  • (4  >F(  II  |T«1»  It) 

»?An*,J 

IFJJ.EO.O)  J«5 
CALL  ANGL 
Ts" 

RPINT  10 

RRTNT  1 1,  IKCR(T) ,1*1 » 19) 

OPtNT  ^ 
fO*  C 

»H'ILT«1  , 

FPL  US»0  . 

1 RrAl*,A,P,H,-p,FMH,VO,’H,JJ 

FC=fhh*i,E-7 
TF|A,EQ.8.)  go  fq  so 
IFfA.GT.C.)  rci  TO  ?0 
A«-  A 

FHMLT*»rULT**ll  ./ri) 

FPLUSsIFLUS/KO 
OPINT  i. 

»hi»lt»i  . 

FPLUS»C  . 

<0*0 

♦0  F<*rHM*l.eE 

!Ff  R.  EO  .r  .)  r,c  T^**  ai 
ST'^R?  A/p 
I*-  *1 

PC*P‘'E0UIM,cO,A  ,q) 

CALL  COKSTAMT,  t9?«A  , 1.  A-fGI 
CALL  HAT^Sfl.H) 

CALL  BlN(>l.fH,AH) 

9r?P«*STNlYC»P’ 

tF<(?I  .€0.0,  ) GO  FO  O? 

oq,oi 

roRp  Cl *9“/F^ 

-^orr  *c  P 

'■R*OTEL?CT<rH,l  ,o,M) 

F»rM 

CALL  C0KSTANfT.lor.A.«.A7f6) 

CALL  MATTSd.H) 
call  pi  a (•!«  ,)«,  SH) 

PX»RM*SINJTC*P) 

FFfPI.EC.O.)  GO  TO  9? 

FI»*(4,,lTi6eE6T:-t/M»1.S»?r»-M*7.i,bfc0?6G€7)/IF6*»<,S?8H.A?.751S*H4 
1M*H**,  ,907T)  ) 

PIsPMFI* 


COPP  CF*P'^/AI 
pH«pM*ccpP€r?**'t» 

FC* ASlN  ISOPTfSI ./PHI  I /P  i 

F.PC»t.€-* 

F*CO»P€C? 

*PLUS*»PIUS*F 
FMULT*rH«'LT*» 

K0*«''*1 

IFTF.lt.I.)  f*1./* 

Tt*9M*  <«tN(  tvr*.Pt)»e)*»’-SI‘lf!TC -.81)*  €)••?)•. 5 
PRINT  ?,A,9,lT0®,H,r«PFC.T0,F,FHH,ER,I0,PN,CCRRECi,RI,CCPRfC?,T€ • 

1 t.RHjTC.HFCf  JJ) 

IFtHOOlI.J)  .t1.  8)  PRINT  4 
GO  TT  1 

9?  PRINT  » 

GO  T)  90 

91  PRIN*  * 

7 roeHAT(7flf,r9.4),lV,RT.t,iv,rE.«i,lV,FS.F,lI,FS.T,?llX,F8.i),lX, 

1FS.?,?MX,F7.1)  •]X,rs.T,l  V«F7.1,tX,Ff.S,lX,FS.?,lX,FS.S,lX,F6.0, 

? IT  ,^9.  3,  tX.AF) 

T foohAT|SX,*A«,*X,»9»,a»,* AT0P»,5X,*H»,SX,*ER«**X,*T8*,4X,*F|HME)* 

t,S»,*FePH?) T,?!, •€••,€ t,« PIN*, SX,*RIN*, FX,*N£AS/*,9X,PR1NP,EX, 

7 •h?AS/*,FX,*.01»X  V<liR  R EDGE  NfN  9F  C*  ,/ , J5X  ,*SPCC*  , 

3  11*,*CALC*.9X,*HEAS*,?»»*RtOAO*,SX,*CALC*,fcX,*M?AS*,9X,*CALC*,1X, 
k •CALC*,fcX,*CAL'**,tt,*TP»OR*,lGX,*»8*,/l 

4 FOOHATl  IM  I 

5 FO'HATC*  E«ROP  9*8*) 

( rORHATl*  EAROP  RlN«e*) 

7 FOOMATl*!*) 

9 F0PHATIA8X,*»HJLT«*,F7,S,SX,*AV£,«* ,77. S,/) 

18  F0OHAT|*Q*) 

11  FooNATllPlAP)) 

99  tTlP 
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1 

9 

• TOT 

N 

t^fC 

ft 

f 

Cti'* 

f 1 *Hf% 
O'; 

*0 

tM 
*'C  05 

» 15/ 

f.otr 

kTH 

4 A a7 

riLC 

.2188 

V .4» 

2 f r ' . 

<•*.* 

VC 

s.orc 

1.4)09 

8.630 

• 3319 

8.90 

.5:9 

123<^.9 

1840.0 

8.45 

iic.c.o 

1.971 

01.6 

1.241 

8.  H 

l.?4l 

4'1  . 

,32- 

* 

i.f:co 

S.C31 

.0316 

8.si 

1.5:9 

i8J*«i 

1849.0 

8.49 

4.i>*l 

iH.ao) 

1 .*«7 

2'. 7. 9 

1.180 

8.74 

t.)70 

4U  . 

\ 1 

S.fOld 

1.0:99 

8.089 

• 8319 

f.19 

i.iro 

t831*9 

1810.1 

8«  •» 

4.4*6 

4‘'0.0)9 

1.494 

302,9 

1.163 

2.73 

1.  US 

4a  • 

.SI, 

U 

1.9'*I9 

1.898 

• 0316 

8.99 

.169 

1?30»4 

1840.3 

8. *8 

8.448 

os.?)! 

1.474 

36.9 

1.190 

1.49 

t.  lie 

171. 

• 54i 

31 

J.9:69 

1.603 

.8310 

8.10 

t.cco 

1834.4 

1841,8 

2.)4 

2.4«9 

136. :0) 

1.159 

111.3 

1.110 

l.*>4 

1*111 

17<  . 

.144 

51 

1.  .'je; 

1.669 

.9316 

8.99 

1.1C2 

180r.4 

1848.1 

»,31 

8.448 

181.030 

1.549 

116.8 

1.161 

1.4) 

1.141 

17*  . 

.*46 

31 

1.1960 

1.190 

.0110 

8.*  8 

.810 

4:o.i 

0*0.) 

8.2* 

2*445 

50.030 

1.691 

46.0 

1.617 

1.81 

1.018 

8C*. 

• 850 

31 

«.5S00 

1.  !:•• 

1 .103 

.8316 

8.90 

l.lfO 

0SO.9 

OIC.f 

8.88 

8.441 

15*.60J 

1.711 

130.1 

1.118 

:.8) 

1.118 

26*?. 

.834 

31 

i.:iC9 

1.033 

• 3116 

2.16 

8.219 

OIC.t 

8.88 

8.44S 

19C,000 

1.101 

108.0 

1.014 

1.14 

1.B14 

1*0. 

.87? 

3 1 

1.9100 

1 .411 

• 1818 

8.10 

.8*9 

4)4.9 

4)6.0 

7.42 

0.087 

1)6.030 

.014 

99.0 

1.013 

1.01 

I.OIC 

4:8  , 

.*‘1  » 

3 1 

V.lrO'JO 

i.  non 

1.0C9 

.1810 

8.11) 

1.463 

413,4 

Ol»i.  4 

2.40 

4,945 

31<J.C0) 

.9)1 

/9).1 

1.16) 

1.1) 

1.  tbi 

4»7, 

• 444 

).l 

1.  9J80 

1.43) 

.1818 

8.f  Q 

8.818 

413.9 

9t9»! 

2. '*8 

1.981 

4*6.660 

.042 

4)3.8 

1.061 

1.91 

1.811 

4'1. 

.5.3 

J1 

i.ec'*: 

».<i)OC 

.►O* 

.1816 

8.93 

.i:o 

It  8J.1 

U40«( 

8.47 

9.114 

*7.413 

.100 

13.8 

.10: 

1.17 

).1»1 

ll**. 

.*10 

3 1 

1.1439 

<|.0?03 

.tif 

•1816 

8.18 

l.OCO 

118*.» 

1191.6 

8.43 

6.813 

18.510 

.298 

163.1 

• 504 

1.11 

1.488 

117. 

• 715 

31 

J.Bl'.O 

%.9190 

.008 

.1866 

8.S8 

1.109 

1171.4 

1144.0 

8. *9 

9.  <54 

ItC.OIQ 

,800 

820.0 

.408 

1.01 

2.07  3 

11). 

.*4? 

31 

i.i9:i 

9.  9398 

1.996 

.1806 

8.10 

.1(0 

ll«4.1 

1194. • 

2.48 

0.  <«:9 

15.000 

• 51) 

4).4 

.917 

1.91 

1.113 

248. 

• 475 

31 

1.9)91 

1.83) 

• 1816 

8.10 

l.rto 

1104, < 

18:4.6 

••*»* 

9.1,*1 

1'1.0)0 

.120 

193.3 

• 901 

1.11 

1.1)1 

24'  . 

.-68 

31 

i.eiio 

1.0109 

1.0)1 

• U%6 

8.10 

1.100 

114».5 

1814. r 

8.  «9 

9.114 

871.&Q8 

.400 

849.3 

.4)5 

1.98 

1.198 

826  . 

.417 

3 t 

l.OS'Q 

t.0)09 

8.093 

• 1810 

8.96 

.8M 

1213.2 

1218. e 

2.41 

16.103 

11. COO 

.0  94 

10.3 

1.18% 

S*10 

1.184 

711. 

.75  7 

31 

1.9931 

t.*990 

8.893 

.1819 

8. 1C 

• 1.8 

1811.2 

18)9.3 

2.4» 

14.519 

8il.6:0 

.87. 

111.0 

1.324 

4.01 

1.324 

87S  . 

.242 

31 

l.J9*e 

1.9983 

1.933 

.1816 

8.18 

.939 

1194.9 

tt09.4 

8.40 

O.OlO 

83.40) 

,471 

29.8 

.15? 

l.*1 

1.244 

82?. 

.41) 

i.s:)9 

•.viei 

1.811 

• 1818 

2.10 

• 340 

Il9(  .1 

1 199. f 

2.4l 

10.983 

29.130 

.954 

30.9 

• HI 

2.31 

1.013 

387. 

31 

1»99)9 

9.3:98 

1.1.9 

.1849 

8.18 

• 8«3 

ti03.>4 

I26>.* 

2.49 

18,0'9 

27,9)9 

.15* 

89.4 

.44,4 

7.71 

1.610 

427. 

.351 

)1 

i.4c:i 

t.  'too 

1.71« 

• 1816 

8.18 

.890 

126*. 0 

18.9.8 

2.*/j 

14.74C 

31.731 

.947 

89.9 

1.101 

3.39 
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41#  . 

31 

1.C939 
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8.013 

• 189C 
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• 819 
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2.49 

14.6^1 
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• 114 
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31 

i.:c9j 
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• 0480 
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• 113 
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843, 
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iri^ti 
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t 
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1811.3 
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1.519 
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i.ce:« 

1.0981 

8.o;o 

• 6318 
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1034.0 

10:9.8 

2. *8 

11*268 

49.909 

.990 

98.i 

1.131 

9.31 

1,131 

912. 

• 101 

Ft 

1.99. 0 

1.  M98 

1 .893 

• 0119 

1.96 

• 801 

|43:.9 

1428.0 

2.18 

4.3*9 

19.7)8 

1.204 

39*1 

1.941 

3.18 

l.f4| 

291. 

.204 

rt 

9.C999 

1.1131 

1.993 

.8319 

1.90 

• ?89 

1041.4 

1980.6 

2.14 

4.119 

19.00I 

1.198 

39.1 

1.888 

3*98 

1.(28 

320. 

• 241 

rt 

9.93;i 

1.1^81 

1 .Ti; 

• 9116 

8.99 

• 242 

ioi;.9 

1438.4 

8.49 

9.984 

99.630 

1.299 

fl.i 

1.098 

4*84 

s.oir 

429. 

• 224 

rt 

|9.93!Q 

0. iiac 

1 .00^ 

.1816 

2.26 

l.ftO 

S9>.« 

203.8 

8.19 

4.000 

34.9)I 

1.439 

J9«l 

1«#80 

• 00 

1.689 

309. 

1.104 

m» 

19.94;9 

0.  1.98 

1.40* 

• •080 

8.89 

1«812 

349.0 

399.9 

8 •t 

2«188 

.f9«80l 

1.989 

J9»f 

l.lOf 

• 98 

1.809 

8«l. 

1.400 

01* 

fo.oaio 

0.4000 

1 .408 

•0310 

8«84 

1.88Q 

393.6 

300.1 

1«331 

14.999 

i.tu 

3i«9 

t.932 

• 19 

1.038 

381. 

8.144 

0 
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PR0GP4M  9INVSH  f OUTPUT • INPUTI 

ON  NOt  IC41.01<S  P4P£R  ON  NICROSTRIP  ANTENNAS  1975 
C*»*TmS  PROSRAH  Pt.OTS  THE  INPUT  RESISTANCE  AT  RESONANCE  VERSES  THE  SUBSTRATE 

c»**thicknfss  chi.  input  data  is  as  eollons. 

C«««F«ERtLOSSrAN,HMlNfHHAX,NH*ABHlN*ABMAX tNARtlL*  VltVE 
C#**F.FR£0  IN  HERTZ,  tR-OIELtCTRIC  CONSTANT  LOSSTAN-LOSS  TANGENT 
C»«*hMTN-NININUN  H TO  BE  PITTED  HMAX-HAXIHUH  H TO  PE  PLOTTED 
C**»NM-NUMB"?  of  HEIGHTS  IPTtNTSI  TO  PLOTTED 

t***  NCTE. CHNAX-hHIN) /INH-1> * CONVENIENT  NUNBER  FOR  THE  H INTERVAL  AND 
C***  (HNAX-HNINI/Zl»  sTONVlfNIENT  NUHBEJ^  FO®  A NICE  GRAPH  SCALE. 

C#*#AP“IN-HtNTNUH  A/B  RATIO  APHAX-M AXIHUH  A/R  RATIO 
C***nA«.NUM9: R OF  RATIOS  (LiNESI  TO  Bt  PLOTTED  ON  THE  GRAPH. 

C***LL-  A LINE  IS  SKIPPED  APTf®  EACH  GROUP  OF  LL  LINES  IS  TYPED 
r»**Vl-MlNrNUN  value  OF  V AXIS  ON  the  GRAPH, 

C***V?-HAX  VALUE  OF  Y AXIS  TN  GRAPH  |RTN  NAX) 

C***  NOTE.  IY>-Y1)/16  * CONVENIENT  INT-RVAL  FO®  G®APH  SCALE.  OR  VlsYZ=e  FO®  AUTO 
C***  ECAlrlNG.  USUALLY  A FULLER  GRAPH  RESULTS  IF  Y1  AND  V?  AFE  OETERHINED  BY 
r***  hand 
c 

R-AL  LAHOA.LAHDAG.L.LOSSTAN.IH 
C0NNnN/9/A,n,F,M,eR 

CO*<HON/r/RS,FC,Pl2,  IM,  L • PI » LA  HO  AC  , L AK- DA  .L  OS  STAN  , ZO.  PA  , PP  . P 
COHHON/F/rFF,OrLTAF .GAIN 
OiMfNSI ON  RPfSl » .X(51» 

«on<i,ji-=i-a/j)»j 

FR-OIUH.E®.  A,9>  - E.R:i^?7165EQ/(  A*SQ®T  ( :®-l»  •(P/HI  ••  .1155» 

l^?.*H*SO®T<P»)l 
IN^I. 

call  ancl 

1 R€  AT*  . ‘P.L^SSTAN,  HHI  N . H^A  X,  NH  , AHH  I‘;  , ABH  AX  .N  A«  ,LL»  Y1  , V Z 

IP(F'^F(qLlNPun  .N-I.OI  GO  TO  R9 
(F  I 

WL  = 1 .1 3/  («"•  SC5T(EC  I » 

NL^l, 

H“IN=MiTN*HL 

MHAX=MMl)f  ml 
WlMCs (HHAX-WHINI/CNM-I » 

ABI^ri(iftMAy-AnMlN»/<NAH-ll 
ATOP=  ABHTN 
CO  2'  Jsl.NAP 
W = MI-IN 
Fr=F*l.-.q 

POINT  ‘^.PF,AT0n,L0SSTAN.'O,HL 

e FO®MflT(//,HY,*FprO(GwT|=*,lPG12,^, JX,*A/n  CATIO  =*,1PGI?.4. 

I TX.*L0S^TAN  =• ,1PG1?. U.Tk,**  o-*,iogi 2.4,3X ••WL  =• .IPGl?. 41 
P5TNT  i. 

U FOOHAT(TX,*A*.ieX,*B*,li,X,*H‘  ,13X,*P  I N^lfT  • , 1 “X  .♦  G A IN  • , 

li:x,*-rr:ci  "NCY»,  7X  .‘T  ANnHlOTM*  .'^X,  -L  A**OAG*  > 

CO  1.  TsI.Nm 

CALL  L-  N'MHfATnp.ZOl 

VO-A/2. 

call  CON^TANC  I,  1 9>  1 , 4’E‘‘» 

call  WATTSd.HI 
CALL  RTN(-l .,H, CH» 

FGsF* I , E-o 

FlXr(w.l31F,f657t-.3/H^i,c«22^M»7.wesC,.‘»66*»7|/(Fr,**(.XZ«i4^2.7515*H4 
t «..o'^7*M*Hn 
p1scm«$IN(vo*p| •♦Z 
RT=FT»FIX 
RO(Tl  =®T 
X(I»sH 

P^TN'f  X.ft.O.H  ,OI,GA1N.-FF,0  LTAF,LA“04G 
3 FORHATlAflPGlfe.c) I 

TFtC-OOd.LLI  .F0.3I  .ASO.IiNH-n  .GT,  (LL/2n>  PC*!  NT  1425 


1425  FOOHATdH  I 

IC  H=H4MlNC 

T1=nh*3,/3.4.5 

TL  OPE  = ( Po  (I  II  -ORClIdNL  /(f  TT»1  I • HINT  > 

POINT  2. glop: 

2 F0ohaT(/.5X,*SL0P£**.pe.,0.*0HHS/NAV  ’ LENGTHS*  I 

TF(J.:').1I  CALL  CALCHl  INH.x,op,r,MMit  ,MHAX,Y1.Y2.12,  ,0,. 

1 tOMKIfoOSTPIP  ANTENNAS, -19, ^MHdSCH' SI, q.QMlNPUT  ^ES. 9.1, 101 
IFfJ.Nf.ll  call  calc**!  (-Nm.X.RC  .-JI 

2C  ATOPaAT^n^ABlNC 

CALL  GotO(j,,T,,.3,,F,ZH,10l 
call  CALOHUC.'’,! 

GO  TO  1 
99  STOP 

FNO 
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lA  a>  rvj 
O'  J 90 
> f fVJ 

to 
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^ ^ r\i  O* 
O'  tc  ^ tst 
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SURROUTIMi  ANGt 

THIS  PCUTINE:  STOP'S  ALL  OP  Thp  VALUSS  FOR  TH=  SIN^  ANO  CCSlhE  FUNCTIONS 
THAT  ARE  USED  THPOUGHT  THE  PROGRAM 
COMMON/A/CO(iAS|,SI(li.5),  ANGINC 
ANGINC=<*. 

OO  10  I=l,l<.5 

ang=anginc*  (i-r t) 

CO( I)=CCS(ANG) 

10  SI( I» =SIN(AKG) 

RETURN 

FNO 


SUBROUTINE  CQMSTBN(MU,SIGHA1 

C THIS  routine  CALCULIT'S  HOST  OF  THE  CONSTANTS  USED  THROUGHT  THE  PROGRAM 

PEAL  K.LAHOAtLAHOAGfML.L.LOSSTAK  ,IM 
COMMON/A/COd-*?)  ,SI(  litE),  ANGINC 
COMMON/0/A,B,F, H,E= 

COMMON/C/PRt’<»P^'*»?P'<»F?A»FCK 
COM MON/C/PCCN 

C0MM0N/E/RS,RC,PI2,IH,L,PI,LAMCAG,LAHDA ,L0SSTAN,70,PA,PP,P 

PI=3. 1415g?e54 

R=63T60 

LAMOAG=2* A*  (4*H/SQRT (cp)) 

LAMOA=1.180Z«5433E10/F 
RS= SORT (PI*F*«U/SIGHfl ) 

RC= RS*A/B 

70=  376.  7303 66 *H/( SORT ( E R ) *3* ( 1 .♦ t . 735* ( EP** (-. 0 724 ))•( H/B 83 6 I ) 

K=PI2/LAM0A 

P=oI?/LAHOAG 

P?=P*P 

L=  70/(F*LAH0AG> 

PCON=ZO/8.* (IM/LAMOA) *•? 

ANGINC=4. 36332313c-? 

PA= A*P* .5 
PSAsP*SIN(PA) 

PC<=K*CCS(PA) 

PB=P*B*  .S 
PSB=P*SIN(Pe) 

PBt<=K*CCS(PP) 

RETURN 

ENT 
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FUNCTION  OIELECT(FN,a,T,iJ) 

this  function  CALCLLflTPS  THE  GIELECTRIC  CONSTANT  REOUIREO  TO  GIVE  THE 
FaEOUfNCV  OBTfllNFO  IN  AN  ACTUAL  TESTED  ANTENNA 
Fh  = "EASUREC  FREOUEKCr  IN  HZ 
A = the  length  a IK  INCHED 
B = THE  NIOTH  0 IN  IKFHES 
H = TH=  SUBSTRATE  TH1''KNESS  IN  INCHES 
F='^  M»  1,  E-9 
BH: .01*  (P/H)**. 11 5E 
fh- (R*H)**2 
ftF= (A  *F) *♦? 

X = T4. 8?6eA?5A*( BH-1) « 

7= ( 4*  FH-AR*FH) 

Y = y*T-278.6l47>»TF*FH 

7 = t»7 
Y=- V/ Z 

0IEL‘CT-Y-SCRT(Y*Y-X*X/Z» 
retu^  n 

END 


SUDFOUTINE  length (AT0B,N) 

this  routine  uses  an  OPTT-IZEO  newton  ITERATION  METHOD  TO  CALCULATE 
7A#  AND  for  a given  A/B  RATIO 

IF  B DOES  NOT  EOUAL  C.O  WHEN  THE  ROUTINE  IS  CALLED.  THEN  B HILL  BE 
USED  AS  THE  FIRST  APPROXIMATION  FOR  B. 

ATOB  = THE  RATIO  A/B 
N = MAXIMUM  NUMBER  OF  ITERATIONS 
COMMON/B/A.B.F, H, ER 

IF(P.eQ,D,»  n=5 .75E9/CF* ATOB* SORT (ER) ) 

K=0 

SER=SQPT(ERI 
FRM=FR-1. 

PSERM=. 61*ERM 
CERM=.0 15  2275*£RM 
T0P  = 5.9!!14?7165E9/F-2*H*SER 
10  01=n 

K=K*1 

OH=B/H 

POT 2= 1, ♦RSFRM’BH**. 1155 
B0T=SQRT(B0T2) 

DBOT=CrRM/(BH** .8845*90T» 

FUNC=TDP/BOT-B*  ATOB 
OFUNC=(-TOP)*OBOT/BOT2-ATOB 
B=B1-FUNC/(0FUNC*1. 334) 

IF(K.GT.N)  GO  TO  23 
IF(ARS(1.-R1/B) .GT.l.E-9)  GO  TO  13 
A=B*ATOB 
RETURN 

20  A=e*ATOB 

15  FORMAK*  ERROR  IN  LENGTH  A =*  , 1PG2  0 . 1 C , * B=*.  1PG2  C . 10  . 

1 • ERROR  IN  B =*,1PGZ0.10» 

PRINT  15,A,B,FUNC 
RETURN 
END 
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SUgPOUT  INE  PAKCUT,  ANGL  £ .OUT,  SC  AlEI 

C TMIS  routine  P900UGFS  A POLAR  PATTERN  ON  THE  CALCOMP  PLOTTER. 

C CUT  and. angle  determine  the  PATTERN  RLANE  THROUGH  THE  ANTENNA. 

C CUTaitY  MEANS  THAT  THE  PLANE  IS  PARALLEL  TO  THE  XY  AXIS  AND. 

C <ANGL“<aTHETA 

C rnp  example,  a plane  going  through  the  X-Y  axis  - ANGLE*90  degrees 
r CUT  a7YX  MEANS  THA T THE  PLANE  PASSES  THROUGH  THE  Z AXIS  ANO 
C <ANGL‘‘<  = PHI 

r FOR  rXAMPLP,  A PLANE  GOING  THROUGH  THE  7-Y  AXIS  - ANGLE^SO  OR 
C 7-X  AXIS  angles  C degrees 

C the  plot  is  scaled  so  that  4C  OB  PLOTTED  IS  ACTUALLY  <SCALE<  D9 • 

C TP  5SCALE<  IS  LESS  THAN  0 WHEN  <PAT<  IS  CALLED.  THEN  THE  VALUE  FOR  SCALE  IS 
r CHANGED  TO  th'  maximum  GAIN  OF  THE  ANTENNA  PATTERN. 

C SCALP  .LT.  0 AUTO  SCALE 

C scalp  .r,T,  0 hOOP  s SCAL*' 

C OuT.L’'.0  PRINT  ANO  GRAPH 

C OUT=l  PRINT  only 

r.  niJTr’  '»CINT  SCALT  ONLY 

r OUT=^  NO  PRT^'TING 

common/ A /CO  (14C  > , SI  (14G ) . ANGINC 

COMPON/PAT/FOXYX ( 77  > ,:pxYYCT’),ETXYX(75» .ETKYY  t 73 > , ETOT X YX f 7S J ,ETO 
t TXVY ( ) 

COhmoM/G/PATk.XV.ZYX 

DIMENSION  ^TMXY  (73»  .jPMXY  <7  3»  .tTvOTXYt  7 7) 
rP(OUT’.GT.?,>  GO  TO  R3 
IF  (CUT. -0 .XY)  PCINT  1 

1 FORMAT(u:x,*X  . Y 3 L a N E*  ,/.1 CX .•phi* ,loX  ,*ephi*, i^x.^etheta 

1*.1?X,*-T0’'AL*,/» 

TP  (CUT.  PO  . 7YX)  P5INT  n*: 

1C'‘  POona  T ( (4?  X, ‘X  - 7 >LfiNE*  ,/ *•  Tm -:T  A*  , 16X  ,♦  ?PHI  •,  lux,  •EThET  A*  ,1  3X  , 

I •-rTOTfiL*,/! 
np  TPsANGL  *73 

IT=IP 
I = C 

XYHAXsC  «l ,r7P| 

no  I"  <1  = 1.  luS, ? 

TP(CUT. -n.tv)  ir=Kl 
IF  (CUT. -O.TYXI  TTs<l 
T=-IM 

UHaUf IT .TP> 

TM=T(TT  » 

lTmxy<:»sPATK*(Um*ST(Tp»*C^(TTM**2 

f'PHXV(:iTpaT<*(U«*CO(TPMTH*‘;l<IT))»«? 

PTnT)IY(T)aPPMXY(I».-THXV(Il 
IFfVYMax.LT.LTOfXYl Tn  XVHiY  = ET1TXY(I» 

1C  CONTTNJ* 

XYOpa  1 ■ ,♦  AL  OG  I?  ( X V'<  AVI 

Tr(*‘r  AL'^.r.T."  ,1  XYl'^aSCALc 

SCAl t aXY^P 

noroNxy  S..C.  -xYpp 

TP  (PUT,  L^  • T .1  PPINT  »,XVOf» 

•»  FOPvaTl*  scale  hC  ot  =•  ,1P'’.2C 

T*** 

TP(PUT.G“.?.»  GP  TO  ?l 

•f  no  ?^  inal,l4S,? 

latM 

ir  (f  PHX  V(  I»  ,NE,  , I ' OHX  Y (11=  l'*  .*  AlOGK  ( - ^^XY  (I  M ♦OBCONXY 

TF  (F  thxyi  II  ,NE. ' . » : thx  Y dial:  .•  ALnr.i  (^ynXY  (in  ♦ORPONXY 

IF  Cf  TOTxy  (T  I .N-  • TOTXYCII  S ic  ,*  ALCGIO  (••  T0TXY(  III  ♦pf'COFXY 

PHla*  .*  (T  -I  » - I*  . 

IF  (r>uT,  L ' . T . > PR  INT  Z,  ^ PmXY  ( I ) , 'T  hxv  ( I ) , : TO  TXY 1 1 1 

9 FfjPHft  T|  U(  lOGl  S.  PM 

IF (fTOTXY (Tl .LT .C  .1  iTOTXYC  Tl aC  • 

TFCFThXYI  II  .LT.C.»-:ThXV  (I|aC. 


rcCtPMXYllI  .LT.C.I  PPHXYCIIaE. 

FPXYXIIlaCailP)  TPHXYI  I) 

PPXY  Y(  Tl  sS  I ( T**  I •£»HXY  ( Tl 
PTXYX(II=C0CIP|  •ETHXYdl 
FTXYY  (I  MSI  (TPI  •pthyy(  1 1 
FTOTX  YX  ( 1 1 am  (ini  • TPT  XY(  II 
FTOTXYYdl  tST  dPI*iTOTXY(ll 

?C  CONTINU: 

IF(OUT.E0.l.»  GO  TO  Zl 

call  CALCHI  I73,pTOTYYX,‘‘TOTXYV,C*?.C.O.D,<J,  .0,  , lOHMl  CRPSTR  IP  ANTE 
lNNAP,-iq,lH  ,1,1H  ,1,1, m 

call  CALCMi c-73 ,epx yx, :pxyy ,*ZI 

CALL  CALPHI (-73,ETXYX, e Txyv,.f I 

call  CALPHKC.O.I 

?1  RETURN 

RNO 
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SUflPOUTINE  RIN(OUT,W,RHI 

THIS  SUBROUTINE  COMPUTES  THE  INPUT  RESISTANCE  AT  RESONANCE  UNADJUSTED 

FOR  the  input  point  - YO 

OUT  .LT.  C.  FOR  NO  PRINTOUT 

OUT  .GT.  C.  FOR  PRINTOUT  OF  RESULTS 

W = WATTS  OF  RAOIATEO  POWER 

PH  = ThR  UNACJUSTEO  input  RESISTANCE  AT  RESONANCE. 
real  LAMOA.LAMOAG.L.LOSSTAN.IM 
COMMON/A/CO(145),SI (145»,ANGINC 
COMMON/B/A, P.F, H.ER 

C0NH0N/E/RS,RC.PI2,  IN,  L , PI , L A MO  AG , L AM OA , L OS  ST AN  ,ZO, PA.PP.P 
COMMON/ F/EFF, CELT AF .GAIN 
COMMON/ G/PATK,XV,ZYX 

NAMFLIST/R£SULT/Z0,L,EFF,W,QR,QC,00,0T,DELTAF,n,  RI  ,RA,RC,RS,F, 

1 LAMDA,LAMDAG,IM,LOSSTAN,A,R,YO,h,ER,GAIN 
OMEGA=PI2*F 
RA=2*W/ (IM* IM) 

QP=  (OMFGA*L*A)/ (2*PA> 

OC=(PI»  70*0 »/(LAMnAG*  RSI 

Q0=1./L0SST AN 

OT=l./( l./OR*l./QC*l./Qn» 

EFF=QT/QR 

DELTAF=F/QT 

RT=RA^2,*RC 

0=  (70*PI/RA»*  {(  (SIN  (PA)  ♦SIN  (PP)  »/  (L  AM0A*PH  **2) 

GAIN=10 ,* ALOGIO  (D*EFF) 

PH=2. •Zn**2/RT 

0 PATK=  rFF*PI* ZO* ( IM/LAM0A)**2/( 2. ’W ) 

IF  (OUT.GT.O.)  PRINT  RESULT 

RETURN 

ENO 


FUNCTION  RIMPd.JJ 

THIS  FUNCTION  GENERATES  thf  PROPER  CORFICIENT  FOR  SUMMING  TERMS  IN  A 
SIMPSONS  RULE  INTEGRATION 
J = TOTAL  number  0'^  TERMS  IN  THE  SUM 
I = TH^  number  I OF  THE  KTH  TERM  3EING  SUMEP 
siMPr2. 

A=I/7. 

11=1/2 

IF(II,E0.A)SIMP=4, 

IF((I,LT.?).OR. (I.GE.J))STMF=1, 

PETU=>N 

ENO 
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I FUNCTION  T(IT) 

THIS  FUNCTION  COMPUTES  A CONVIENIENT  INTERMFOIATE  NUMBER 
IT  - refers  TO  A PARTICULAR  ANGLE  FOR  THETA 
REAL  K 

COMMON/ A/CO (145) »SI (145) t ANGINC 
COMMON/B/A,B,F,  H,ER 
COMMON/ C/P2, K,PSB,P8K, PSA, PCK 
BKC=K*B*CO(IT)* .5 
T3=PSB*  COS(BKC) 

T4=PRK*SIN(BKC) *C0( IT) 

T8=P2-<K»CO (IT) )**2 
T=  (T3-T4)/T8 
RETURN 
END 


FUNCTION  U(IT,IP) 

THIS  FUNCTION  COMPUTES  A CONVI^NIENT  INTERMEDIATE  NUMBER 
IT  - REFERS  TO  A PARTICULAR  ANGLE  FOR  THETA 
IP  - REFERS  TO  A PARTICULAR  ANGLE  FOR  PHI 
REAL  K 

COMMON/ A/CO (145), SI (14 5), ANGINC 
COMMON/R/A, B,F, H,£R 
COMMON/C/P2,K,PSB,PBK,PSA,PCi< 

AK=K*A*SI(IT)*SI(IP)* .5 
U2=PSA*COS( AK) 

U3=PCK*SIN(AK)*SI  (IT)*  SKIP) 

U5=P2-(K*SI (IT) *SI(IP) )**2 
U=(UZ-U3)  /U5 
RETURN 
END 
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SUBROUTINE  MATTSINT.M* 

THIS  SUBROUTINE  CONPUTfS  THE  RADIATED  POWER  BY  USING  SIHPSONS  RULE  TO 
INTEGRATE  WITH  PHI,  AND  THEN  USING  FILONXS  HETHOO  TO  INTEGRATE  WITH 
THETA. 

NT  = I FOR  half  sphere  INTEGRATION  I AS  KALOI  SUGGESTS) 

NT  = 2 FOR  FULL  SPHERE  INTEGRATION 
W = WATTS  OF  RADIATED  POWER 
COHNON/A/CO (IASI .SI (IAS) .ANGINC 
COHMON/O/PCON 

OIHENSION  Y (381  ,G(38), AX(38) 

IF((NT.GT,2l.OR.(NT.LT,l)IGO  TO  AO 
IF(CO(BO)  .LT.  .S)  CALL  ANGL 
XANG=2.*ANGINC 
Nl  = 73 
M2=1AS 
W=0. 

N=0 
30  N=N«-1 

J=C 

DO  10  IT=H1,M2,2 
J=J41 

AX(J)=ANGINC*  (IT-73) 

TH=T(IT) 

1=0 

G(U)=0. 

00  20  IP=  37,109,2 
I=I*1 

UH=U(IT,IP) 

20  G(J)=G( J) ♦( (UH*CO(IPV»TH*SI (IT) ) * • 2» ( UH* SI ( IP) • CO (I T) )••2)•SIHP (I, 

1 371 

10  G(J)=G( J)*XANG/3, 

W=W*ABS(XFIL(AX ,G, Y ,X A NG, 37 , 1 . , 1 ) • pCON ) 

Hl  = l 
H2=73 

if(n,lt.nt)  go  to  30 
RETURN 

AO  PRINT  1,NT 

NT  = C 

NT=1/NT»1/NT*1 

RETURN 

1 FORHAT(*  NT=*,  120,*  ERROR  NT  IS  AfJ  INTEGER*) 

END 
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FUNCTION  XFIL(TtF,V,rH,NT,X, Jl 

THIS  OOUTINF  COHPUTSS  THE  INTEGRAL  OF  SIN  OR  COS  FUNCTIONS  OF  THE  FORM 
GIX) -INTEGRAL!  F (T» ‘SIN ! X» T I ) OT  FROM  T1  TO  TN 
USING  FILONXS  method. 

T = ONE  DIMENSIONAL  ARRAY  <0IM  GP'ATFP  THAN  NTI  CONTAINING  VALUES  OF  T 
F = ONE  dimensional  ARRAY  JOIM  GR^^ATER  THAN  NTI  CONTAINING  VALUES  OF 
F(TI 

Y = ONE  dimensional  ARRAY  (DIM  GR‘’ATEP  THAN  NTI  USED  BY  TH  ROUTINE 
TH  r SPACING  BETWEEN  POINTS  IN  RADIANS 
NT  r number  of  POINTS 
X r THE  parameter  in  XSIN<X*TIX 
J = 1 FOR  SIN  FORM 
J = ? FOR  COS  FORM 

dimension  T (II ,F(1)  ,Y( II .alpha (91 .BETA (91 .GAMMA (91 
COMMON/A/CO (1491 .SI (145I.ANGINC 

DATA  (AL^HAdl.  1 = 1 . 91  / 0 . . - . 03  0 C 02  30  CO  0 1 724.  . C 0 0 0 3 00  00 143  307  . 

1 -.0  CO  00  00  093P683E  ..OC  3 0 00  4E9P  41  71..-.  CO  00  171 02239324*.  00  0 4232804232 

2 8.  >.0  9 8X493  0E.Xl,a?:6..  3 4444444  444444  4/.  (BE  T A ( 1 1 . I:  1, 91 /- . 0 00  0 00  0 00  0 

3 32325..  29 00 00 00 2395272  . -.00 03 00  13X8  2.344. .cocoa  5923 00592.-. 00 01 795 

4 7351  290  T..  0 0 3*^2  73  35  85  3 6 7.-.  C 7 8 0 9253  8 . 9 523  8 ..  1 33  3333  3333  3 333  , 
5.66E566656666667/.(GAMMA(II .I=1.9l/c..-.c00c00000000180..00000000a 
6 042834. -.  00  or  90  CD  77  083  41..  uOC  001 '■a  208  4335.-.  COO  0881  83421 517..  00  476 
71904761925.-. 133333333333333. 1. 33333333333333/ 

M=  M0D(NT.2I 
IFlM.Ea."!  GO  TC  66 
THrTA=  X»TH 
A=B=G=C . 

IF(THETA.Gc..75I  go  TO  3 
THETA’S  THETA**2 
B=  .OOa000000000349*THETA2 
DO  4 Isl,« 

As  (A.ALRHA (II  I ’THETAa 
Rs  (B*BFTA( II l•TH£TA^ 

4 G=  (G’GAMMAdll  •THETA2 

As  (AtALPHA (911 ‘THE TA*THETA2 
3 H=  P»BETA(91 
G=  G+ gamma (91 

I71=X»T(11*22. 91 83118  ♦73. 

I7NsX*T(NTI*22. 9183118  ♦73. 

IF(J.EQ.ll  GO  TP  6 

7 no  8 K»1.NT 
IZ=X*T(K)*22. 91 83118  ♦73. 

8 Y(Kls  F(KI»CO(I7l 

XFs (F (NTI*S1(1ZNI-F (II  ‘SI (I71))*A 
GO  TO  9 
6 DO  5 K=1«NT 

I2=X*T(KI*22. 91 83118  ♦73. 

5 Y(KI=  F(KI»SId7l 

XFa  (F  (NTI  •CO(I2NI-F  (II  ’COdZl  II  • (-AI 

9 SUMS  (Y (ll^Y(NTII/2. 

Nils  NT-1 

DO  10  Ks3,nT1,2 

10  SUMS  SUM^Y(KI 
XFs  XFtB’SUM 
SUMS  0. 

DO  11  8=?, NT. 2 

11  SUMS  SUM^Y(KI 
XFs  XFtG*SUH 
XFIL*  XF*TH 
RFTIIPN 

66  XFILs  3. 

PRINT  100. NT 

100  FORMATdHO.lTHNUMBER  DF  POI  NT  S=  . 1 10 . 3 6H  MUST  BE  ODD  XFiL  SET  EOUA 
IL  TO  Z-ROI 


RETURN 

END 
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